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SUMMARY 


A shallow spherical dome subjected to lateral pressure is a structure 
for which the deformation departs appreciably from the linear theory at 
relatively small values of the deflection amplitude. It is also one for 
which the buckling process is characterized by a rapid decrease in the 
equilibrium load once the buckling load has been surpassed. For struc- 
tures having this type of buckling characteristics the question arises 
as to whether the proper buckling criterion to apply is the classical 
criterion, which considers equilibrium with respect to infinitesimal 
displacements, or the finite-displacement "energy criterion" proposed by 
Tsien. 

In this paper the problem of the finite displacement and buckling 
of a shallow spherical dome is investigated both theoretically and exper- 
imentally. In the theoretical approach the nonlinear equations are con- 
verted into a sequence of linear equations by expanding all of the vari- 
ables in powers of the center deflection and then equating the coefficients 
of equal powers. The basic parameter for the shallow dome A is propor- 
tional to the ratio of the central height of the dome h to its thick- 
ness t. For small values of this ratio the expansions converge rapidly 

and enough terms are computed to determine the buckling load according 
to the classical criterion. For higher values of h/t, convergence 
deteriorates rapidly and it was not possible to determine the buckling 
load with the number of terms which were computed. However even for 
these higher values of h/t the deflection shapes are determined for 
deflection amplitudes below the amplitude at which buckling occurs. These 
deflection shapes are characterized by their rapid change as h/t 
increases and by the fact that, over most of the range of h/t studied, 

the maximum deflection does not occur at the center of the dome. 

Experimental results seem to indicate that the classical criterion 
of buckling is applicable to very shallow spherical domes for which the 
theoretical calculation was made. A transition to energy criterion for 
higher domes is also indicated. 
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INTRODUCTION 


The development of the theory of bending of thin-walled spherical 
shells has a long record. A survey of the problem can be found in refer- 
ences 1 and 2. The fundamental equations are developed by Hans Reissner 1 
(1912) who shows that, for a thin-walled spherical dome that is not shal- 
low, the membrane stresses in the shell maintain equilibrium with the 
external pressure, while the bending of the shell has relatively little 
effect except near the edge of the shell where the shell adjusts itself 
quickly to the prescribed boundary conditions. Bending in the shell is 
therefore essentially an "edge effect" or "boundary layer" phenomenon. 
Asymptotic solutions of the bending problem have been obtained by 
Blumenthal (1912), Havers (1935), Jacobsen (1937), and others 1 on the 
basis that the parameter (R/t) 2 is very large, where R is the radius 
of the spherical shell and t, its wall thickness. Both symmetrical 
and nonsymmetrical loading and edge conditions have been treated, 
including the case of a dome supported on columns. 

The asymptotic solutions are, however, not valid for shallow spheri- 
cal shells, 2 for which the effect of edge conditions is no longer limited 
to a thin layer near the edge and the interaction of bending and membrane 
stresses is strong. In 19*4-6 Eric Reissner (ref. 3) developed the gov- 
erning equations for shallow spherical shells on the explicit assumption 
that the ratio h/a is so small that (h/a) 2 is negligible in compari- 
son with h/a, h being the height of the dome and a its radius (see 
fig. 1). A few special cases are solved in reference 3« 

Reissner' s solutions are based on linearized equations. Since the 
effect of bending on the membrane stresses is strong in the case of a 
shallow dome, one naturally asks the question: To what extent is the 

process of linearization valid? Expressed in terms of the ratio of the 
vertical deflection at the center of the dome to the wall thickness w Q /t^ 
the question is: How soon does the solution deviate from linearity as 

w Q ^t increases? 

To answer this question the nonlinear problem is treated in the 
present paper. The particular problem of a shallow spherical shell with 
a clamped edge carrying a uniform pressure is chosen so that a convenient 
experimental comparison can be made. It is shown that the nonlinear 
character depends upon a parameter A which is defined as 



^ee Timoshenko's book, reference 1, for references to original papers. 

2 By shallow is meant a spherical segment for which the ratio of the 
height to the base radius is small, say, less than l/8. 
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where p is the Poisson's ratio and t, a, and R are as previously- 
defined (see fig. 1). The range of w 0 /t in which the linear solution 
is valid is small indeed. For example, at A = 4, the equilibrium pres- 
sure given by the linear solution is, respectively, 9, 23, and 50 percent 
too high when w Q /t is 0.1, 0.25, and 0.5. 

Consider now the problem of buckling of thin-walled spherical shells. 
For a complete sphere under uniform pressure, the classical solution, on 
the basis of linearized equations, is obtained by Zoelly (1915), Schwerin 
(1922), and Van der Neut (1932). (See ref. 2, p. 491.) The buckling 
stress a cr is given by 



where q cr is the critical value of external pressure. This stress has 

the same magnitude as the critical stress for an axially compressed 
cylindrical shell of radius R and of thickness t. It is relatively 
high in comparison with experimental results. The corresponding buckling 
mode predicted by the theory is also at variance with laboratory experi- 
ence. To reconcile the differences between theory and experiment Von 
Karman and Tsien in 1939 (ref. 4) introduced a new concept into the 
theory of elasticity: the "lower buckling load." They discovered that 

for values of pressure q considerably below that given by equation (2) 
quite different stable states of equilibrium exist, which could be 
revealed only by abandoning the classical linearization of the problem. 
The minimum of such values of q is the lower buckling load q^. If 
q exceeds q^, the chances are great that buckling will occur. In 
reference 4 the lower buckling load is computed (subject to a number of 
simplifying assumptions) with respect to a special class of buckling 
modes. Friedrichs in reference 5 avoids some of the arbitrary assump- 
tions by applying asymptotic integration in the manner of a boundary- 
layer theory. Application of Friedrichs' equations, however, yields no 
minimum buckling load, and it is pointed out (ref. 6) that the minimum 
obtained in reference 4 is due to the special form of displacements 
assumed in that investigation. 

The final "energy criterion" of buckling is formulated by Tsien in 
reference 6. It is stated that under average laboratory and actual serv- 
ice conditions the most probable equilibrium state is the state with the 
lowest possible energy level. In other words it is assumed that there 
are disturbances of sufficient magnitude so that the transitions from 
higher energy levels to lower energy levels are always possible. Two 
conditions must be satisfied in defining the "possible energy levels": 

(l) the corresponding external forces and internal stresses must be in 
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equilibrium; (2) the geometric restraint and loading conditions, if any, 
must be satisfied. Tsien points out that these necessary conditions for 
possible energy levels are not checked in references 4 and 5* When the 
check is applied (ref. 6), a lower buckling load is obtained for spheri- 
cal shells on the basis of Friedrichs * equations. The agreement with 
experiments is good. 

It appears that these arguments apply equally well to spherical 
domes and to the complete sphere. Therefore, the first theoretical ques- 
tion to be settled is whether the "classical criterion" of buckling or 
the "energy criterion" should be used in calculating the critical buckling 
load. The classical buckling criterion is based on the assumption that a 
given state of equilibrium of a shell becomes unstable when there are 
equilibrium positions infinitesimally near to that state of equilibrium 
under the same external load. Thus in applying the classical criterion 
an equilibrium state is compared with its neighboring equilibrium states 
and the incipient buckling is revealed by a negative slope of the load- 
deflection curve, that is, when an increase of deflection corresponds to 
a decrease in the corresponding applied load. The important contrast 
between the classical criterion and the energy criterion is that in the 
former only a continuous load-deflection process is considered, while in 
the latter a jump to the state of lower energy level is permitted even 
though the intervening states involve higher energy levels . The linear- 
ization of the governing equations, ordinarily made purely for mathemati- 
cal simplicity, should not be regarded as a part of the classical 
criterion. 

Although the energy criterion seems plausible, nevertheless it can 
be verified only by comparison with experiments. The energy criterion 
necessarily yields a buckling load which is never greater than that given 
by the classical criterion. If there is a wide difference between the 
two buckling loads the problem becomes simply to choose the criterion 
that gives closer agreement with the experiments. 

For shallow spherical domes the buckling load calculated on the 
basis of the classical criterion, but without linearizing the governing 
equations, is known only in very few cases. In the comparison with 
experiments presented in figure 20 of refererce 6, the curve labeled 
"classical theory" is really the one given by equation (2), which is 
applicable to a complete sphere and is calculated from linearized equa- 
tions. When the nonlinear equations applicable to a shallow spherical 
dome are used the buckling load is lower than that given by equation (2). 
For example, when A = 4 the calculation of the present report gives a 
buckling load which is about one-half that given by equation (2). Thus 
the wide difference between the classical theory and experiments exhibited 
in the figure cited above may be entirely caused by an improper mathemati- 
cal process. 
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To clarify the argument further, consider the case of a flat arch, 
as a two-dimensional analog of the spherical dome. For such an arch two 
buckling modes are possible. If the arch rise is high, it buckles in the 
mode shown in figure 2: the center line of the arch remains essentially 

" inext ensional." If the arch rise is small, it may buckle downward with 
a sudden reversal of curvature, as shown in figure 3: a phenomenon some- 

times described as "oil-canning" or "durchschlag." The axial compressive 
strain plays a dominant role in the latter case and linearization of the 
governing equations is not permissible. A detailed study made in refer- 
ence 7 shows that in practice the classical criterion agrees better with 
experiments, except for very low arches (arches whose rise is of the order 
of the wall thickness) for which the energy hump tends to vanish and the 
gap between the two criteria tends to be closed. 


For shallow spherical domes the prevailing buckling mode is of the 
oil-canning type, in which the membrane stress plays an important part, 
and is basically a nonlinear phenomenon. 


There exists only one paper on the oil-canning of shallow spherical 
domes based on the classical criterion. This is Biezeno's work (ref. 8) 
which treats a shallow dome whose edge is free to expand so that the mem- 
brane stress in the radial direction vanishes on the edgej and the dome 
is subjected to a concentrated load acting at the center. The following 
equations (which are equivalent to those of the present paper) are 
obtained: 


r2 




v 0 + r‘ 




1 ~ jj r% 2 _ Pr (l - h 1 2 ) /l dA 

2 T 2 n Et \R dry 


( 3 ) 



where P is the central load, \|r = — is the slope of the deflection 
_ . dr 

surface in a meridional section, w being the radial displacement normal 
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to the original spherical shell, and v Q is the component of displace- 
ment normal to the axis of symmetry, that is (see fig. l), 


Other symbols are defined in the list which follows this section. 

Biezeno makes the following simplifying assumptions to obtain a solution: 

(1) that the term on the right-hand side of equation (3) may be neglected; 

(2) that in equations ( 3 ), (5) , and on the right-hand side of equation (4) 
the slope of the radial displacements \|/ can be written as 


where C-^ and C 2 are two undetermined constants. Equation (4) is then 

solved with proper boundary conditions. Let the solution be denoted by 
^ 2 > which, of course, is different in form from equation ( 7 ). Biezeno 

then determines the constants and C 2 by requiring that typ and ^2 

yield the same values of vertical displacement at the edge of the plate 
(r = a) and at the center (r = 0). The load-deflection curve can then be 
calculated from equation ( 7 ) and the buckling load determined. 

The influence of Biezeno' s simplifying assumptions on the buckling 
load is not easy to assess; and there exist no experimental results to 
compare with the theory. 

The case considered in the present paper is that of a shell clamped 
at the edge and subjected to uniform lateral pressure. The equations of 
equilibrium (equivalent to eqs. ( 3 ), (4), and ( 5 )) are solved as pertur- 
bation series expressed in powers of the parameter w Q /t, that is, the 
ratio of the deflection on the axis of symmetry and the wall thickness 
of the shell. The load-deflection curve so determined is used to obtain 
the buckling load. 

Relatively few assumptions are made in the present calculation. 
Unfortunately the perturbation series seems to deteriorate rapidly for 
large values of A, so the result is satisfactory only for A of order 5 
or smaller. In this range of A the buckling loads computed on the 
basis of the classical criterion agree quite well with experiments. 

On the other hand the calculation of the buckling load on the basis 
of Tsien's energy criterion also offers considerable difficulty. If the 
formulas of reference 6 are extended to cover the shallow shells studied 
in the present paper it is found that the so-called "lower buckling load" 


v Q = u cos 9 - w sin 0 


( 6 ) 


* = *3. = C 1 | + c 2 | log f 


(7) 
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has an equal or higher value than that given in equation (2) when A < 10. 
This unreasonable result is obtained because the energy expressions and 
the mode shape assumed are not sufficiently accurate. It is not clear 
how to improve the results. Theoretical deflection curves derived from 
the bending theory do not permit a very simple representation. In any 
case, the refinement of Tsien's calculation would have been a major 
endeavor. For the same reason the calculation of the buckling load on 
the basis of classical criterion using the Rayleigh-Ritz method is not 
pursued. Therefore the most convenient theoretical determination of the 
critical buckling load remains an open question. 

In order to do justice to either the classical criterion or the 
energy criterion, further theoretical study must be made for A in the 
range, say, from 5 to 15 . A study of all available experimental data on 
the subject seems to show that the classical criterion of buckling holds 
for very shallow spherical domes, while a transition to energy criterion 
takes place at some intermediate values of A of order 6. 

One more point should be mentioned before the presentation of the 
main analysis. In reference 5> Friedrichs suggests that it may be possi- 
ble that a boundary layer occurs at the edge of a certain segment , the 
width of which in its turn shrinks to zero with the thickness of the 
shell. This suggestion seems plausible because as the shell becomes 
thinner and thinner the bending of the shell becomes less and less impor- 
tant . In the limit t — > 0 the deflected surface must be an "applicable" 
surface of the original.-' In the upper part of figure 4 the shell repre- 
sented by the dotted line is applicable to that represented by the solid 
line; in other words, a deformation of the solid line into the dotted 
line involves no strain energy due to the membrane stresses. To account 
for the small but finite bending energy of the shell the deflection sur- 
face may take the form represented by the lower part of figure 4. A 
boundary layer may be developed at the segment angle a. This conjecture, 
however, turns out to be improbable for a shell subjected to uniform 
external pressure; since it can be shown that the segment angle a tends 
to zero at a higher order in t (the shell thickness) than does the 
boundary- layer thickness. Therefore the boundary layer can be developed 
only at a pole a = 0 which is the case presented in reference 5. 

The investigation presented in the present paper was conducted at 
the California Institute of Technology under the sponsorship and with 
the financial assistance of the National Advisory Committee for Aeronautics. 


5 Two surfaces are called "applicable" to each other in differential 
geometry if one can deform into the other by continuous bending without 
stretching or tearing the surface. 
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SYMBOLS 




integration constants 


a base radius of shell 

a ,b n , c n coefficients in power-series expansion of F 2 in terms 

of Ax; see equation (A5) 


D = 





E 

E n 

fn 

h 

K 


coefficients in power-series expansion of F 2 in terms 
of Ax; see equation (A7) 

Young's modulus 

functions of f n and w n ; see equation (28) 
coefficient of expansion for S r in powers of W Q 
central height of shell above base plane 
constant; see equation (Al) 


k = ^12(l - p 2 ) 


“t 

N r 


Pn 


radial bending moment per unit length 
circumferential bending moment per unit length 
radial membrane force per unit length 
circumferential membrane force per unit length 


E 


5 A 

; (t) q 


coefficient of expansion for P in powers of W Q 


Q shear force per unit length perpendicular to middle sur- 

face of shell 

q pressure on surface of shell; positive when directed 

downward 


t 
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R 


initial radius of curvature of shell 


r 


horizontal distance from axis of symmetry of shell 




t 

u 


W = w/t 

W o = w oA 


thickness of shell 

radial displacement of middle surface of shell measured 
tangential to initial surface and positive in outward 
direction 


vertical displacement of middle surface of shell meas- 
ured perpendicular to initial surface and positive in 
downward direction 


w n coefficient of ejqpansion for W in powers of W Q 

w Q = w(r = 0) 


x = r/a 


z 


o 


a 

a n>Pn->7'n 

P 

6 


initial distance of point on middle surface of shell 
above base plane 

segment angle of a possible deflected surface 
functions of 9 n and 9 n j see equations (35) 

semi- included angle of shell 
finite-difference interval 


e r radial strain 

6 t 


circumferential strain 


10 


NACA TN 5212 


T] = X^/4 
0 = sin”l(r/R) 

0 n particular integral of equations ( 2 6 ) for f. 

* 2 - # - - >* 2 ) I 

li. Poisson's ratio 

p = bei A ber'A + (l - ber A) ber'A 

cp circumferential position angle 

cp n integrals of w n j see equations ( 29 ) 

Subscript : 

cr critical 


THEORETICAL ANALYSIS 
Derivation of Equations 


Consider the spherical shell segment of radius R, base diameter 2a, 
height h, and constant thickness t shown in figure 1. The initial 
position of a point in the central surface is given by the cylindrical 
coordinates r, cp, and z 0 , where r is the radial distance from the 
center, measured parallel to the base, cp is the circumferential angle, 
and z Q is the vertical distance, measured upward, from the base plane. 

It is assumed that h/a is small enough that 


z Q = h + \/r 2 - r 2 - R ~ h - | y 

Asp ___ r 
dr R 


( 8 ) 


The deformation of the middle surface is assumed to be radially 
symmetric and is therefore specified by u, measured tangential to the 
middle surface in the outward radial direction, and w, measured perpen- 
dicular to the middle surface in the downward direction. The deflections 
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are considered to be finite, but small enough so that 
neglected with respect to unity. 



can be 


Since, under these conditions, the magnitudes of vectors tangential 
to the middle plane are 'equal to their components parallel to the base 
plane, the equations for forces and moments in the middle plane are 
identical with those for a flat plate. That is, 


d_ 

dr 


^rM r ) - Mj. - rQ = 0 


( 9 ) 



( 10 ) 


where N r and N-^ are, respectively, the radial and circumferential 
membrane stresses, M r and are the corresponding bending moments, 

and Q, is the shear stress in the direction perpendicular to the deformed 
middle surface. Vertical equilibrium of a central cylindrical section of 
radius r (fig. 5) requires that 

r9 ^ + " r ( 3? 2 • 3 ?) (11) 

where q is the applied pressure. Substituting equation (ll) into equa- 
tion ( 9 ) and using the approximation equations (8) result in 






rq dr + rN- 


Jl + ^L)=o 

\R dr ) 


( 12 ) 


The bending moments are expressed in terms of the deflections using 
the strain-deflection relations 





( 13 ) 
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where e r and are the longitudinal strains of the middle surface 

in the radial and circumferential directions, respectively. Then 



M r 



if ^ t 

r dr 




(Ik) 




-D 


/l dw 
y r dr 


+ 



where 





Using these expressions for 


M r and M_^ , equation (12) becomes 


Dr 


jL 1 

dr r dr \ dr J 



(15) 


This is the first basic equation. Now from equations (13) 


u 

r 




+ 


w 

R 


so that 


du 

dr 


-L JL/W - prN„ + Et 
Et dr( t r 



(16) 


These values are substituted in the first 
a second relation between N^. and N r : 


of equations (l4) to obtain 
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N r = 


Et 


— 2 

A- JL^n. - |irN r + Et SSL') - 2 + + JL/n, - jjlN ) 

Et dr \ ^ r R/R 2\dr/ Et\ t V 


( 17 ) 


Combining this equation with equation (lO), the second basic equation 
is obtained: 


2 

— i — (r 2 N r '') + i Etf — ^ + Et £ ^2 = 0 
dr r dr \ r / 2 \dr/ R dr 


( 18 ) 


Knowing w and N r , can be obtained from equation (10) . In 

the problem to be studied q is a constant so that equation ( 15 ) becomes 


D A. 
dr 


r dr \ dr / r dr / 2 


qr 


( 19 ) 


These equations are transformed into nondime ns ional form by the use 
of the following variables: 


x = £ 


W = 2 


k 2 = 12(l - 1 1 2 ) 

,2 


S-r = 


N 


3 

Et 


St = 

Et'' 


P = 


l 2 / \^ 

1 


E 


: (t) 




,2 _ ka 2 
A “ Rt 


( 20 ) 
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The parameter A^ can also be expressed as (see fig. 1) 

A^ = k ^ sin^p « k 


or 


2 h sin E _ k h (1 + cos p, 
t 1 - cos p a 


2k i 


( 21 ) 


O 

Thus for the assumed range of (3, A is proportional to the ratio 
of the central height of the dome to its thickness and can therefore be 
interpreted as representing the ratio of the compression stiffness to 
the bending stiffness. 

Upon substituting these new variables, equations (l8), (19), and (10) 
become 


— i — (x 2 S r ) 

dx x dx' 



2l ^ 

k dx 


0 


( 22 ) 


I A. 

x dx 


1 _d_/ x dW \ 
x dx\ dx / 


k A 


+ 


k m s 


X dx/ r 


6 p 


(23) 



(2k) 


With A = 0 (R = oo) these are Harman' s equations for the finite 
deflection of a flat plate, expressed in polar coordinates. Their deri- 
vation is exactly analogous to Chien's derivation of the equations for 
the finite deflection of a flat circular plate (ref. 9)* 


Expansion in Terms of W Q (=w 0 /t) 

As in Chien's paper (ref. 9) the procedure used for solving equa- 
tions (22), (23), and (24) is to consider the center deflection 
ratio W(0) = W Q as a parameter and to expand all of the variables in 
powers of W Q . Thus 
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P = p i W o + p 2 W o 2 + p 3 w o 5 + • • • 

2 3 

W = WjCxJWq + w 2 (x)W q + Wj(x)w 0 + . . . > 
S r = fi(x)W 0 + f 2 (x)W 0 2 + f 3 (x)W 0 3 + . . . 


(25) 


These expansions are valid for small enough values of the deflection 
ratio W Q , but their exact range of convergence is unknown. For the 
case of the flat circular plate Chien obtained good convergence for values 
of W Q as high as 4. 

Substitution of these series in equations (22) and (23) and the 
equating of equal powers of W Q result in a sequence of pairs of simul- 
taneous equations for f and w n . Each of these pairs of equations 
can then be combined to obtain an equation for f n alone 

1 _d_ 

x dx 


_1 _d_ ^ 
x 3 dx ' 


X' 


5 _d_/l df n ) 
dx\x dx J 


^ xK _ 6.2 

>•+ Af n ---Ap n 


n 


(26a) 


or 


dV 


dx 


ip + i 

4 X 


d5f. 


dx" 


x^ dx c 


n 


x^ dx 


S + A 4 f 


n 


- ! v 2 - F „ 


(26b) 


plus an equation for w n in terms of f n 

= + <27) 
where Eh is a constant and F n and cp n are the following functions: 
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F-l = 0 


F3 = kA 



1 d i/ dw i\ 2 


1 dx x dx 

X dx 2^dx j 


(f 0 ^i + f, *sUj4- 

d / dv l dwg\ 

\ 2 dx ^ dx / x dxj* c 

Lx\dx dx / 


v _ kX 2 1 L ^3 + f ^2 + f £^1 ) + 

F k ' kA xl fl ST +t 2d5T +1 3dx / + 


1 _d_ 

x dx 


1 _d_ 
< x dx 


dw^ 

dx 


dv 3 if^V 

dx 2\dx / 


(28) 


cp x = 0 

cp 2 = 


^3 = 


94 = 


r 

l/ dw l\ 2 

dx 

J 0 

2\dx J 

X 

r 

dwj_ dw 2 

dx 

J 0 

dx dx 

X 

r 

l/ dw 2 \ 2 

( 

+ . 

Jo 

2\dx / 

< 


dwj 

dx” 


dx 


(29) 


Upon making the substitution br\ = x 2 equation (26a) becomes 




A 2r jf n 


" ' an x2pn 



(30) 


which can be recognized as the equation for the lateral deflection of a 
linearly tapered beam on an elastic support whose spring constant is a 
linear function of position along the span. This interpretation is use- 
ful in the numerical work which follows. 
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The solution of the homogeneous part of equation (26a) or (26b), that 
is, the complementary function, is 

f n = ^A n her* Ax + B n bei 1 Ax + C n ker f Ax + D n kei*Ax^ 


where 


ber 1 z = — ber z 
dz 

and the ber and bei functions are defined in terms of J Q , the zero- 
order Bessel function of the first kind, by 

J G ^zi^/ 2 ) = ber z + i bei z 

with an analogous relation between ker z, kei z, and 

Boundary Conditions 

The boundary conditions for a clamped-edge shell subjected to a 
radially symmetric distributed load are 


at x = 0, 


at x = 1, 


— = 0, S r is finite 
dx* 


W = p. = U = 0 

dx 


(31) 


To satisfy the first two conditions it is necessary that 


Cn - Dn - 0 

In terms of the expansion coefficients, equations (25), the remaining 
boundary conditions become 
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w 


n 


= 0 




dwn 

dx 


= 0 


x = l 


(1 


n)f n + x 


d£n 

dx 



(32a) 


Because of the nature of the expansion there is the add.it ional con- 
dition that 


w n(°) 



n = 1 
n ^ 2 


(32b) 


The constant E n in equation ( 27 ) can be eliminated by combining 
equations (32a) and (32b) so that the boundary condition on w n becomes 


w n (0) 



n = 1 
n > 2 


(32c) 


Let 0 n be the particular solution corresponding to F n> on the 
right-hand side of equation (26a) ; then the complete solution for f n 
becomes 


f n = ^( A n ber,7oc + B n bei'Xx) + 0 n - | (33) 

while 

w n = ^( A n be i ^ x " B n ber Ax ) " ^2^ — ( x ^ 9 n) "" ^2 + F n 

Substitution of these values into boundary conditions (32) and their 
solution for A n , B n , and p n result in 

A n = ber'A + p n (l - ber A)j ( 3 ^a) 
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B n = bei '* " Pn bei >) 


Pn = 


kh c 


6(1 - n)p 


■|(l + p) J^(ber' A) 2 + (bei'A)^J 


A(ber’A bei A - bei'A ber A) 

r r 


+ 


(i + p) 0 1 - ber A) ber ’A - bel A bel'xj + 

1 ^ Pn^ + ?n 


where 


p = bei A ber'A + (l - ber A) ber'A 


„ _ 1 
^ A 


1 A( x 2 e ) 

x dx\ n / 


X=1 v2 

+ — <P n (l) 
„ p-\ n' • 

x=0 


r* 

— - 

-N 

( — 

£1* 

; s( x2e ") 

> 

J 


x=l 


y _ kA 2 (e + 1 d O 

■ ~^\ a ~ trj x=1 


First-Order Solution 

The particular solution 0q of the first-order equation is 
and the equations for Aq, Bq, and pq reduce to 


(34b) 

(34c) 

(35a) 

(55b) 

(55c) 

(55*) 

zero 
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A 

A ber'A 
kp 



A(ber'A bei A - bei'A ber A) r 


The values of A-^, and are given in table 1, while the 

values of and f-j_ are given in tables 2 and 3 and. are plotted in 

figures 6 and 7, respectively. This first-order solution is identical 
with the linear solution previously found by Reissner (ref. 3). 

For the higher order equations no solution was found in terms of 
known functions and so it was necessary to resort to power-series expan- 
sion and numerical methods. 

Power-Series Solutions 

Judging from the work of Chien (ref. 9) it was felt that calculation 
of the first two terms , p^ and P 2 > of the expansion for the pres- 
sure P would permit at least an approximate determination of the 
buckling load. Therefore a power-series solution for was obtained 

even though it was realized that the succeeding solutions could not be 
obtained by this method because of the involved form of the functions F n . 
The procedure and formulas used are shown in appendix A. Since the 
expansions are all in terms of Ax it was recessary to restrict the 
calculations to values of A = 8. The values of P 2 obtained are shown 

in table 4 and plotted in figure 8. These values ax-e negative for small 
values of A, but become positive at A = 6.5 and are rapidly increasing 
at A = 8. Since buckling can occur only when some of the p n 's are nega- 
tive it was clearly necessary to obtain the higher order terms of p n . 

These additional values of p n were obtained numerically. 

Numerical Solutions 

The differential equations (26) contain the unknown parameter p n 
and also have the unwieldy boundary condition that S r is finite at 

x = 0. As a result a complete numerical solution would be very difficult. 
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Therefore only the particular solution is determined numerically, using 
arbitrary boundary conditions. The required boundary conditions are 
then satisfied using the known solution of the homogeneous equation. 

In terms of the finite-difference approximations 


f (x) = ^-[f(x + 5) - f(x - 6)] 
f"(x) = -^[f(x + 6) - 2f(x) + f(x - 5)] 

f'"(x) = -i-[f(x + 26) - 2f(x + 6) + 2f(x - 6) - f(x - 26)1 

? fiP l J 


(37) 


flV ( x ) = ^[f(x + 2 &) “ + 5 ) + 6 f(x) - Af (x - 6 ) + f(x - 26 )J 

where 5 is the difference interval; equation ( 26 a) I with the constant 

3 Pn ^ 2 \ ' 

term - f- omitted becomes 

8 k / 

(\ + ^ 9 n( x ) “ ( -j- + -I i-x + — ^)e (x + 6 ) + 

\6^ & 2 x2 ) n \ b k & 3 X 6 2 x 2 26^?/ n 

( ' 55 + + 552 + si?) 0 n< x - 5 > + + gf;) 9 + + 


^ ' is) 6 "'* ' 25) = l6Fn(x) 


(38) 


The desirable boundary conditions for 0 n are the ones which give 
a smooth solution. For the tapered-beam analogy of equation (30) the 
obvious boundary conditions meeting these requirements are those repre- 
senting unsupported ends. When these boundary conditions are transformed 
in terms of the variable x they become 
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at x = 0, 


at x = 1, 


d0. 

dx" 


n _ 


da 


n 


d-'e 


n 


dx^ 


d J B r 


dx^ dx5 


= 0 


= 0 


(39) 


A first attempt to solve the finite-difference equations (39) by- 
relaxation was unsuccessful because of slow and erratic convergence. 

Instead, Crout's method of solving simultaneous equations (ref. 10) was 
used to determine the values of 9 n at 11 points (& = 0.10) at once. 

This could be done rapidly, but unfortunately 11 points were not enough 
to determine accurately the end values and derivatives which were 
required. Instead of decreasing the spacing to 0.05 throughout, it was 
decided to add two end sections from 0 to 0.3 and from 0.7 to 1.0 with 
0.05 spacing. The solutions in these end sections were joined with the 
original solution at the 0.30 and 0.70 stations where the function and 
its first derivative were matched. Since the higher derivatives were 
small at the junction points (especially for A = 4 and 7) this method 
was adequate, but did cause some trouble when higher derivatives were 
required for the succeeding calculations. 

Calculations of 82 were made for A = 4, 7, 10, and 13, while 
for A = 4 and 7 the calculations were continued to determine 8j 

and 0^. As A increases convergence of the series for P deteriorates 

rapidly and the function 0 n has increasingly large oscillations. It 

was decided therefore not to continue the calculations for A = 10 and 13. 
The values of p n obtained are shown in table 4, while the values of w n 

and f n are shown in tables 2 and 3 an( l plotted in figures 9 an( l 10. 

For A = 4 the convergence of the series for p n was very satis- 
factory, the contribution of the fourth-order term still being small at 
the critical buckling load. In figure 11 is shown a plot of load P 
versus the center deflection ratio W Q and in figure 12 are shown the 
deflection modes for several values of the center deflection ratio W Q . 
These deflection modes have their maximum at the center, and as W Q grows 
they become increasingly peaked toward the center. 

For A = 7, however, the convergence is poor and the coefficients p n 
all being positive, no buckling can be determined using just four terms. 
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The convergence is good enough to determine the deflection shapes W 
for small values of W D and these are plotted in figure 12. These 
deflection shapes give an explanation of why p 2 is positive (which 
implies increasing stiffness with respect to the center deflection as 
the load increases) since they show that the maximum deflection is no 
longer at the center and that with increasing load the center deflection 
becomes a progressively small portion of the maximum deflection. This 
characteristic is corroborated by the experimental measurements. 

The deflection modes for A = 10, which are also shown in figure 12, 
exhibit the same characteristic, but with the position of the maximum 
deflection moved Outward toward the edge. However, since these curves 
are calculated using only two terms of the expansion for W, these curves 
should not be taken as indicating accurately what happens at the larger 
values of W Q . 

The rapid change which must occur at buckling from a shape in which 
the maximum deflection occurs near the edge to one in which the max imum • 
deflection occurs at the center is probably also an explanation for the 
poor convergence. Since the experimental results show that at A = 10 
the maximum deflection is again at the center it may be that for these 
higher values of A the convergence is actually improved. However, to 
obtain accurate values of 0 n for these higher values of A it would 

be necessary to start with a smaller finite-difference interval than was 
used here. 

Since the influence of the p^ and pjj. terms at the buckling load 
for A = 4 was small, it was felt that for A < 5 an adequate approxi- 
mation to the buckling load could be obtained using just the first two 
terms p^ and p 2 . The critical conditions occur when "Pt - =0 so 
that for ° 


P = Pl W o + P2 W o 2 


the critical conditions are 


W r 


-Pi 


cr 2p 2 
2 


cr 


"Pi 

kp 2 




> 


(^o) 


2k 
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This value of P cr is plotted in figure 13 where it is compared 

wrch the experimental results. The minimum value of A for which P cr 

exists is that for which the critical deflection equals the initial 

height of the dome, that is, for W = h/t. 

u cr 

EXPERIMENTAL PROGRAM 
Equipment 


An experimental program was carried out on a series of shallow domes 
having a base diameter of 8 inches, nominal radii of curvature of 20 
and 30 inches, and nominal thicknesses varying from 0.032 to 0.102 inch. 

The edges of the specimens were held between two rings which were bolted 
to a circular plate (figs. 14 and 15) thus providing a rigid built-in edge 
support and a closed pressure chamber. A separate set of clamping rings 
was used for each of the two radii of curvature. The specimens were 
subjected to a uniform normal load using both oil and air pressure; the 
oil provided an approximation to a constant volume characteristic during 
buckling while the air provided a constant pressure characteristic. 

The specimens were made by spinning from flat sheet. After unsuccess- 
ful attempts to heat-treat aluminum spinnings, magnesium alloy QQ-M-44 was 
selected because of its favorable ratio of yield stress to Young's modulus 
compared with other non-heat -treated metals. Magnesium also has the 
advantage that since it is spun while hot most of the residual stresses 
are eliminated. This is evidenced by the small separation when a radial 
cut is made in a magnesium spinning. Because of the difficulty of 
spinning such shallow shells the preliminary specimens were very dis- 
appointing, but by a combination of spinning on concave and convex molds 
the quality was greatly improved. Unfortunately it is still not so good 
as would be desired. 

Pressure measurements were made using a Bourdon tube for pressure 
over 20 psi and a mercury manometer for pressures under 20 psi. Excep- 
tions were two of the early specimens having low buckling loads which 
were tested using the Bourdon gage. This gage of course gives a closer 
approximation to a constant volume characteristic than does the manometer. 

Deflection measurements were taken with a 0.001- inch- scale dial gage 
riding on a channel beam fastened at its ends to a. circular ring which 
rotated in a groove cut in the upper clamping ring. Readings were made 
to the nearest 0.0005 inch. Traverses were made on two or more diameters 
to determine the initial shape of the shell and were repeated at intervals 
during the loading. Intermediate measurements were also made of the cen- 
ter deflection. Because of the variations of the specimens from a true 
spherical form the question arose as to what should be taken as the 
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radius R from which the parameter A was calculated. It was decided 
to assume that the central rise h would determine the radius since A 
can be simply expressed in terms of h (eqs. (21)) and because experi- 
ence with the buckling of shallow arches showed that for arches having 
the same central height small symmetrical variations in shape have only 
a small effect on the buckling load. In figure 1 6 the variations from 
the assumed radii are shown for typical examples of each of the six com- 
binations of the two nominal radii, 20 and 50 inches, and each nominal 
initial sheet thickness 0.032, 0.054, and 0.102. It is seen that the 
variations increased markedly with the thinness of the sheet and the 
flatness of the dome. 


Oil Tests 

The oil pressure tests were made first, and two or more tests of 
each combination of thickness and radius were made. The early preliminary 
tests made on aluminum samples all showed a very distinct unsymmetrical 
buckling mode. This is believed due to the high residual stresses 
resulting from the spinning operation since the majority of the magnesium 
specimens buckled symmetrically. In the cases in which unsymmetrical 
buckle did occur in the magnesium specimens the mode was not of the over- 
all unsymmetrical form such as the unsymmetrical mode of vibration of a 
flat circular plate. Rather it appeared that the buckles themselves were 
inherently symmetrical but were displaced from a central position on the 
shell, probably because of initial asymmetries of the shell. 

The unsymmetrical buckling occurred only in the range of A between 
6.0 and 8.6 and was associated with a prebuckling deflection mode in 
which the displacement at about half the radius from the center was 
greater than that at the center. 

In figure 16 are shown the deflection curves of the specimens. 

There is a distinct change in the deflection modes as A increases. 

^ near 4 the deflection is peaked at the center and decreases 
steadily toward the edge. As A increases, the peak gradually flattens 
out, until at A = 5-^5 the maximum deflection no longer occurs at the 
center. Instead, at large- deflections there are two peaks symmetrically 
placed at about a half radius from the center. With a further increase 
in A the peaks move outward until finally when A = 8.98 a third peak 
appears in the center . This gradually becomes the predominant peak. 

These trends agree very well with the theoretical deflection curves for 
A = 4, 7, and 10 shown in figure 12. 

In figures 11 and 17 are plotted the curves of pressure versus cen- 
ter deflection of the specimens. For low values of A (A< 5), the 
specimens buckled in a continuous manner. As more oil was pumped into 
the chamber the pressure increased more slowly, reached a maximum, and 
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then decreased. But for A > 5> the process was discontinuous. Usually 
there would be a slight movement of the shell without the addition of oil 
followed by a sudden jump to a lower pressure and a greater displacement. 
There was no regular trend in the ratio of the pressure after buckling 
^P cr j to the buckling pressure P cr as A increased and also sur- 
prisingly no significant difference in this ratio between the tests made 
using the Bourdon tube and those using the manometer. 


Air Pressure Tests 

For the air pressure tests an accumulator tank was connected to the 
air line close to the testing fixture so that the buckling process was 
practically a constant pressure process. Buckling occurred very suddenly 
and with a sharp report. The final buckled shapes were symmetrical with 
deflections very much larger than those of the oil tests. Deflection 
traverses were made during loading, but it was inadvisable to make them 
at loads approaching the expected, buckling load. Two examples of these 
deflection traverses are shown in figures 16 (g) and 16 (h). 

The specimens remained in their buckled position after the pressure 
was released. An approximate determination of the pressure required to 
unbuckle them was made by unbolting the clamping rings, inverting the 
rings with the specimens still placed between them, and then bolting the 
inverted assembly to the base plate. The pressures required to unbuckle 
the specimens were considerable and are included in table 5- 


DISCUSSION 


The physical parameters and buckling loads of all the specimens are 
shown in table 5* In figure 18 are plotted the buckling loads as a func- 
tion of A. The oil pressure tests are shown with black dots, while the 
air tests are shown with open circles. For the oil tests the points at 
the lower ends of the dashed lines indicate the value to which the pres- 
sure jumped during the buckling process, wh/ili a wing on the left of a 
lower circle indicates an unsymmetrical buckling mode. When plotted on 
logarithmic paper the results tended to follow two intersecting lines. 

In figure 18 the corresponding power-law curves are shown. 

In figure 13 the experimental buckling loads are compared with the 
theoretical loads calculated using two terms of the series for A < 5> 
and the one point calculated using four terms for A = 4. In figure 11 
the corresponding curves of load versus center deflection are also com- 
pared. Although the experimental results are low compared with the 
theory, the difference (approximately 15 percent at A = 4) is not great 
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considering the variations of the initial shapes from a true spherical 
surface. Part of the difference can be attributed to yielding which 
occurred at the higher loads , especially for the specimens having values 
of A near 5* It is felt that the results are close enough to corrobo- 
rate the theory proposed and establish the applicability of the classical 
criterion for buckling for the lower range of A (A < 5). 

From figure 18 it is clear that the type of loading, air or oil, has 
little if any effect on the buckling load. This is expected to be the 
case if the classical buckling criterion holds, but would be rather 
unexpected if Tsien ! s "energy criterion" applied. For in Tsien f s theory 
when the buckled and unbuckled energy levels are compared the loss in 
potential energy of the load during buckling must be included. Since 
the loss in potential energy is a maximum when the pressure remains con- 
stant, the buckling in a constant -pres sure system should occur at a lower 
load than for any system in which the pressure decreases during buckling. 
Rather large differences in the buckling loads between a "rigid testing 
machine" (approximated by oil loading) and a "dead-weight loading" 
(approximated by air pressure loading) are predicted by Tsien in refer- 
ence 6, according to the energy criterion, for complete spheres and also 
for spherical domes of fairly large values of A. As remarked before, 
however, the calculations as given in reference 6 are probably inappli- 
cable to domes as shallow as those tested here, and a more accurate energy 
expression and deflection mode should be used. Although it is impertinent 
to reject the energy criterion on the basis of the disappearance of the 
difference in buckling loads between the rigid testing machine and dead- 
weight loading, the fact that the buckling loads given by the classical 
criterion are reasonably close to the experimental values for A < 5 
seems to indicate that, at such small values of A, a refinement of the 
energy-criterion calculations is unwarranted. 

The buckling stress for spherical shells is usually expressed in 
the form 


a 


cr 




2t qcr ) 


(41) 


where K is a numerical constant. The classical theory for the buckling 
of complete spheres gives a value of K of order 0.606 with p = 0.3 
(see eq. (2)). Tsien' s theory in reference 6 gives a K of order about 
half of that given by the classical theory. The calculation based on 
the classical criterion given in this paper, for very shallow domes, also 
gives a value of K ranging from 0.2 to 0.4. Both reference 6 and the 
present paper predict K as a function of the parameter A. A composite 
picture showing the results of all available theoretical and experimental 
data is given in figure 19, where K is plotted against the parameter A. 
It is to be remarked again that all the curves labeled "Theory" should 


28 


NACA TN 3212 


be read with caution. On the extreme left, the theoretical result 
according to equation (40) does not converge for values of A > 5* and 
extrapolation to larger A values is dangerous. On the other hand, 
the curve showing theoretical results of reference 6 probably is good 
only toward the right-hand end; the left-hand end of that curve probably 
should be more carefully computed. In other words, to do justice to 
either the classical criterion of buckling or the energy criterion of 
buckling, both theoretical curves shown in figure 19 should be recalcu- 
lated to higher accuracy. The curve showing the "classical linear 
theory" applies only to complete spheres and is shown here for reference. 

Even though not proved rigorously, a trend of shifting from the 
classical criterion to the energy criterion at A of order 6 seems to 
be indicated by the experimental results. This transitional value of 
A is probably low on account of the imperfections of the test specimens. 

A comparison of the above results with the case of low arches as 
given in reference 7 reveals a fundamental difference, with respect to 
the stability criteria, between the buckling under lateral forces of a 
spherical dome and an arch. In the arch case the classical criterion 
holds for larger values of A (i.e., higher arches), and a transition 
to the energy criterion occurs as the arch becomes very low. In the 
dome case the classical criterion holds for smaller values of A (i.e., 
very shallow domes), and a transition to the energy criterion occurs as 
the dome becomes higher. 


SUMMARY OF RESUITS 


In this paper an attack upon the problem of the finite deflection 
of a shallow spherical shell has been made. The results may be sum- 
marized as follows: 

1. The theoretical approach has been to transform the nonlinear 
equations into a sequence of linear equations by expanding all the 
unknown functions in powers of the nondimens ional center deflection W Q 
and equating coefficients of equal powers of tf 0 . The initial linear 
equation can be solved exactly in terms of the ber and bei functions, 
but the succeeding equations have had to be solved either by power series 
or numerical methods. For small values of the parameter A the resulting 
series converges rapidly enough so that a determination of the buckling 
load can be made using only four terms of the series. For higher values 
of A the convergence deteriorates rapidly, so that for A greater than 
5 no determination of the buckling load can be made. However, for deflec- 
tions smaller than the critical buckling deflection, the deflection modes 
can be determined for a much wider range of A. These deflection modes 
change rapidly with A and for values of A near 7 have the surprising 
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characteristic that the maximum deflection occurs approximately halfway 
between the center and the edge of the shell. 

2. For small values of A the results of the experimental program 
agreed substantially with those of the theoretical analysis. The buckling 
load near A = 4 was only about 15 percent below the theoretical value 
while the trend of the buckling loads as A increased was approximately 
the same as predicted by the theory. The deflection modes also showed 
the same characteristics as predicted by the theory. The experimental 
buckling mode was inherently symmetrical as assumed in the theory; the 

few exceptions can be attributed to large initial asymmetries in the 
specimens . 

3* Tests were made with both air and oil pressure, which approached 
the extremes of constant pressure and constant volume buckling character- 
istics, respectively. The buckling loads obtained by the two methods 

showed no significant difference; thus a feature of the "energy criterion" 
did not appear. 

4. Experimental results seem to indicate that the classical criterion 
of buckling is applicable to very shallow spherical domes and that a tran- 
sition to energy criterion occurs for higher domes. 


California Institute of Technology, 

Pasadena, Calif., August 13, 1953. 
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APPENDIX A 


INFINITE-SERIES EXPANSION 


After substituting the expressions for f n and w n , equation ( 55 ), 
into the second of equations (28) for F2, it becomes 



k 2 /a 2 

32x 2 \d x^ 


1 _d_ 
x dx 



bei A x - Bq ber' Ax 


(Al) 


The particular integral of equation (26a) corresponding to the first 
term of equation (Al) is 


0 


2 



bei 1 Ax - Bq 


ber 1 Ax 


) 


(A 2 ) 


The particular integral for the remaining terms of F2 (all quad- 
ratic in ber' and bei') is obtained by expanding in series. The series 
expansions for ber' Ax and bei' Ax are 


ber 'Ax 


bei 'Ax 


■ (&*? 


\ 4 n-l 

£<->' 

n=l (2n - 1)'. (2n)I 


^4n+l 


n=0 (2n)'.(2n + 1)1 


(A 5 ) 


of 


In terms of the above series the expansion of the quadratic terms 
becomes 
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p , i = _ k 2 A 2 

64 


j^lj^n + b n + 2 ( n + x )( 2n + 1 ) c n +l] (^) 


( A l 2 - B l 2 ) C n + n ( 2n + 1 )( A l 2 a n - B! 2 b n ) (^) 


v 4n-2 


(A4) 


wae re 


a n = (-) 


n 


m 


=0 (2n - 2 m)i( 2 n - 2 m + 1 ) • ( 2 m) I ( 2 m + 1 ) ' 


b n = (") 


n 


m=l (2n - 2m + l)'.(2n - 2m + 2)'. (2m - 1)1 (2m)'. 


► (A5) 


c n = (-) 


n 


m: 


=1 (2n - 2 m)l( 2 n - m + 1 ) 1 (m - 1 ) 1 ml 


The corresponding particular integral of equation (26a) is 


0»i- 


(A6) 


where 


~ " a n-l “ ^n-1 “ ^n-lj 




n > 1 


£>n ” 


(2n - l) 2 [(2n - l) 2 - l] 
c n-l + S n -l] 




l)(2n - l)a n _ x + 


n ^ 2 


(A7) 


h„ = 


(2n - l) 2 |(2n - l) 2 - l] 
n-1 " b n-l| 


i( n - l)(2n - l)b n-1 + 


c 

g Q = h o = 0 


n > 2 
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The values of d , gj_, and hj_ are arbitrary with respect to the 

recursion formula of the differential equation, but when equation (a 6) is 
combined with equation (A2) it is required that 


d 


o 



h 


1 



(A8) 


in order that 0^ = 02 * + 02 * 1 have the proper limiting value as 

X — > 0. 

The coefficient g-^ is completely arbitrary and for convenience 
was taken as equaling unity. 


The infinite-series expansion for cp^ is 


Then ©2 and its derivatives and cp^ are substituted in equa- 
tions (35) to determine and p n which are in turn substituted into 

the boundary conditions, equations (34), to obtain A 2 , B 2 , and P 2 * 



cp 2 


-r 

n=0 


A l 2a n " B l 2b h| 


4n + 2 

l 2 / 


4n+2 


a(r) 


4n 
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TABLE 1 


VALUES OF A-j_, Bj_, AND Pl 


A 

A 1 

B l 

Pl 

0.5 

-4.875 

155.9 

5.333 

• 7 

-5.481 

56.78 

5.361 

1.0 

-2.454 

19.39 

5.441 

1.5 

- 1 . 6 l 4 

5.614 

5.877 

2.0 

-1.194 

2.200 

7.061 

3.0 

-.7298 

.4093 

14.35 

4.0 

-.4275 

-.06698 

35-97 

5-0 

-.1855 

-.2078 

89.55 

6.0 

-.005114 

-.1892 

206.2 

7-0 

.07882 

-.09905 

426.1 

8.0 

.07857 

-.01571 

784.7 

9.0 

.4214 

-.02332 

1,311 

10.0 

.009895 

.02615 

2,043 

11.0 

-.005544 

.01498 

3,031 

12.0 

- . 007745 

. 004468 

4,343 

15.0 

-.004895 

-.0009003 

6,048 

14.0 

-.001756 

-.002132 

8,219 

15.0 

-.00005555 

- . 001506 

10,932 

j 


TABLE 2 


VALUES OF w n FOR A = 4, 7, 10, AND 13 




A = 

= 4 



A = 

7 


A = 

10 

II 

X 

w i 

W 2 

w 3 

w 4 

W 1 

w 2 

w 3 

w 4 

W 1 

W2 

wi 

0 

1.0000 

0 

0 

6 

1.000 

0 

0 

0 

1.000 

0 

1.0000 

• 05 

.9965 

-.0021 

-.0018 

-.0008 

1.0011 

.0080 

.020 

.052 

1.0002 

-.0014 

.9999 

.10 

.9859 

-.0082 

-.0072 

-.0032 

1.0044 

.032 

• 079 

.206 

1.0010 

-.0053 

.9995 

.15 

.9683 

-.0181 

-.016 

-.0068 

1.0093 

.070 

•173 

.453 

1.0025 

-.0091 

.9988 

.20 

.9436 

-.0312 

-.027 

-.011 

1.0152 

.121 

• 300 

.782 

1.0053 

-.011 

.9979 

• 25 

.9116 

-.0470 

-.039 

-.015 

1.0211 

.182 

.449 

1.173 

1.0099 

-.018 

.9970 

• 50 

.8724 

-.065 

-.051 

-.019 

1.0254 

.250 

.614 

1.604 

1.0168 

.017 

.9964 

.40 

.7725 

-.101 

-.073 

-.023 

1.0214 

.386 

.938 

2.435 

1.0381 

.112 

• 9993 

• 50 

.6453 

-.130 

-.081 

-.019 

•9825 

.477 

l.l4l 

2.952 

1.0625 

.300 

1.0148 

.60 

.4958 

-.141 

-.074 

-.0092 

.8832 

.480 

1.129 

2.903 

1.0608 

• 517 

1.0438 

.70 

.5538 

-.126 

-.053 

.0006 

.7016 

.366 

.863 

2.187 

.9708 

.614 

1.0437 

• 75 

• 2533 

-.110 

-.039 

.0034 

.5787 

.289 

.684 

1.777 

.8672 

.581 

• 9955 

.80 

.1769 

-.087 

-.025 

.0037 

.4391 

.191 

.473 

1.242 

• 7149 

.463 

.8851 

.85 

.1085 

-.059 

-.013 

.0031 

.2920 

.106 

.280 

.740 

.5169 

.306 

.6945 

.90 

.0525 

-.032 

-.0048 

.0018 

.1528 

.040 

.127 

.348 

.2938 

.139 

.4290 

• 95 

.0143 

-.0087 

-.0011 

.0052 

.0448 

.0087 

.032 

.089 

.0931 

.034 

.2960 

1.00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


vn 
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TABLE 3 


VALUES OF f n FOR A = 4, 7, 10, AND 13 




A = 

= 4 



A : 

= 7 


II 

H 

O 

roy 
1 — 1 

II 

r< 

X 

f i 

f 2 

f 3 

f 4 

*1 

*2 

f 3 

f 4 

f i 

f l 

0 

-4.244 

1.353 

0.180 

0.0788 

-16.246 

-2.038 

-7.04 

-17.56 

-37.21 

-65.43 

.05 

-4.240 

1.354 

.180 

.0781 

-16.250 

-2.069 

-7.11 

-17.77 

-37.21 

-65.43 

.10 

-4.227 

1.358 

.183 

.0760 

-16.262 

-2.159 

-7.34 

-18.38 

-37.22 

-65.42 

.15 

-4.206 

1.364 

.186 

.0724 

-16.282 

-2.305 

-7.71 

-19.38 

-37.23 

-65.41 

.20 

-4.176 

1.372 

.191 

.0675 

-16.306 

-2.504 

-8.21 

-20.72 

-37.25 

-65.40 

.25 

-4.137 

1.381 

.196 

.0612 

-16.334 

-2.749 

-8.82 

-22.36 

-37.27 

-65.39 

.30 

-4.089 

1.391 

.200 

• 0538 

-16.360 

-3.030 

-9.52 

-24.23 

-37.32 

-65.37 

.40 

-3.968 

1.409 

.207 

.0361 

-16.392 

-3.670 

-11.08 

-28.39 

-37.45 

-65.37 

.50 

-3.811 

1.418 

.207 

.0169 

-16.353 

-4.293 

-12.54 

-32.25 

-37-64 

-65.44 

.60 

-3.623 

1.410 

.198 

-.0002 

-16.178 

-4.760 

-13.58 

-34.92 

-37.82 

-65.67 

• 70 

-3-409 

1.377 

.184 

-.0127 

-15.795 

-4.-952 

-13.92 

-35.80 

-37.76 

-65.95 

• 75 

-3.294 

1.350 

.176 

-.0163 

-15.507 

-4.929 

-13.83 

-35.57 

-37.54 

-65.96 

.80 

-5.177 

1.317 

.167 

-.0179 

-15.153 

-4.837 

-13.58 

-34.94 

-37.12 

-65.73 

.85 

-3.060 

1.278 

.159 

-.0182 

-14.737 

-4.697 

-13.22 

-34.02 

-36.46 

-65 . 09 

.90 

-2.945 

1.236 

.152 

-.0180 

-14.271 

-4.530 

-12.79 

-32.91 

-35.56 

-63.92 

• 95 

-2.837 

1.193 

.146 

-.0175 

-13.780 

-4.360 

-12.33 

-31.75 

-37-46 

-62 . 18 

1.00 

-2.737 

1.153 

.l4l 

-.0170 

-13.304 

-4.206 

-II.90 

-30.65 

-33.30 

-60 . 12 
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TABLE 4 


VALUES OF p n (n > 2) CALCULATED 


A 

P 2 

p 3 

P4 

Power 

series 

Numerical 

1.0 

-1.120 




1.5 

- 2.726 




2.0 

-4.919 




3.0 

-12.39 




4.0 

-26.47 

- 26.5 

3-77 

1.08 

5-0 

-45.2 




6.0 

-38.1 




7-0 

53-7 

56.6 

314 

801 

8.0 

219 




10.0 


392 



13.0 


537 
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TABLE 5 


EXPERIMENTAL DATA 
[„ = 0.52, E = 6.5 X lO^ psij 
(a) Hydraulic pressure tests 


Specimen 

in. 

h, 

in. 

A 

Qcr > 
psi 

P cr 

Method of pressure 
measurement 

(a) 

( P cr ) 2 

1 

0.101 

0.251 

4.04 

36.2 

12.3 

B 


2 

.099 

.253 

4.08 

32.5 

11.8 

B 


5 

.101 

.256 

4.16 

35.8 

12.1 

B 


k 

.100 

.365 

4.80 

60.5 

21.0 

B 


5 

.101 

.376 

4.94 

72.5 

24.8 

B 


6 

.055 

.24-0 

5.45 

14.0 

61.3 

M 

25-3 

7 

.055 

.251 

5-57 

12.1 

54.3 

M 

25-3 

8 

.052 

.297 

6.08 

15.2 

71.3 

B 

48.3 

9 

.055 

.380 

6.75 

31 

122.4 

B 

60.0 

10 

.051 

.4-10 

7.22 

27.5 

147 

B 

85.1 

ll 

.051 

.4-22 

7.40 

25.2 

136 

B 

60.9 

12 

.031 

.303 

8.04 

4.2 

165 

B 

106 

15 

.032 

.361 

8.59 

6.02 

201 

M 

124 

14 

.031 

.353 

8.69 

4.59 

185 

M 

137 

15 

.033 

• 394 

8.82 

7-33 

213 

M 

ill 

16 

.033 

.4-10 

8.98 

8.96 

255 

M 

179 

17 

.029 

.4-44- 

10.1 

6.7 

354 

B 

201 


(b) Air pressure tests 


Specimen 

■t, 

in. 

h, 

in. 

A 

c 4:r> 

psi 

p cr 

Pressure required 
to unbuckle specimen, 
psi 

18 

0.101 

0.382 

4.98 

73.5 

24.9 

16.5 

19 

.101 

.426 

5.26 

99.5 

33-8 

35 

20 

.055 

.265 

5.62 

16.85 

65.2 

— 

21 

.054 

.413 

7.10 

33-6 

143 

12.6 

22 

.032 

.347 

8.45 

5.67 

190 

1.8 

23 

• 033 

.399 

8.91 

11.75 

350 

3.3 


a B, Bourdon gagej M, manometer. 
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Figure 1.- Notation. 





Figure 2.- Buckling mode for a high arch. 



Figure 3«- Possible buckling 
mode of low arch. 
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o 
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Figure 4.- "Applicable” and approx- 
imately "applicable" surfaces. 
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Figure 5«- Equilibrium of a central cylindrical section. 



Figure 6.- Variation of linear deflection mode w x with A. 
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Figure 7 .- Variation of linear membrane stress f 1 with A. 
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Figure 8.- Variation of p-^ and P 2 with X. 




(*) A = 7. 


w n for A = 4 and 7. 


Ul 


Figure 9.- First four terms of 
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ON 




(a) X = 4 . (b) X = 7 - 

Figure 10.- First four terms of f n for X = 4 and 7 . 
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Figure 11.- Comparison of experimental pressure - center-deflection 
curves with theoretical results. 
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Figure 12.- Deflection modes for A = 4, 7, and 10. For A = 4 and 7, 
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Figure 13.- Comparison of experimental buckling loads with theoretical 

values . 
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Figure 14.- Testing jig. 
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IN. FROM CENTER 



DEFLECTION FROM INITIAL POSITION 

VARIATION OF INITIAL POSITION 

FROM A CIRCULAR ARC THROUGH CENTER 



IN. FROM CENTER 


(a) Specimen 1; A = 4.04; 

h = 0.251 inch; P cr = 12.3; 

oil loading. 


(b) Specimen 5; A = 4.94; 

h = 0.376 inch; P cr = 24.8; 

oil loading. 


Figure 16.- Experimental deflection shapes. 
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(c) Specimen 6 ; X = 5.45; 

h = 0.240 inch; P cr = 6 l.J; 

oil loading. 


Figure 16. 


BUCKLED DEFLECTION, IN. 


DEFLECTION FROM INITIAL POSITION 
'VARIATION OF INITIAL POSITION FROM 
.A CIRCULAR ARC THROUGH CENTER 



(d) Spe cimen il; A = 7.40; 

h = 0.422 inch; P cr = 135 . 9 ; 
oil loading. 


Continued. 


STSf VOW 



4 2 0 2 4 

IN. FROM CENTER 


DEFLECTION FROM INITIAL POSITION 



VARIATION OF INITIAL POSITION FROM 
A CIRCULAR ARC THROUGH CENTER 



IN. FROM CENTER 


ui 


(e) Specimen 12; A = 8.59; 

h = O. 56 I inch; P cr = 200.6; 

oil loading. 


(f) Specimen 16 ; A = 8 . 98 ; 

h = 0.L10 inch; P cr = 25^.5; 

oil loading. 


Figure l 6 .- Continued. 
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BUCKLED DEFLECTION, IN. 


DEFLECTION l-ROM INITIAL POSITION . 

(VARIATION OF INITIAL POSITION FROM 
[A CIRCULAR ARC THROUGH CENTER 



(g) Specimen 19 ; A = 5-26; 
h = 0.426 inch; P cr = 33 . 8 ; 

air pressure loading. 


(h) Specimen 21; A = 7.10; 
h = 0.413 inch; P Qr = l43; 

air pressure loading. 


Figure l6.- Concluded. 
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Figure 17.- Experimental pressure - center-deflection curves. 
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Figure 18.- Experimental results. 
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Figure 19«- Comparison of theory and experiment. K = 
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